# -*- coding: utf-8 -*-
# @Create On    : 2020/4/24 16:49
# @Author  : Bao Linfeng
# @Function    : 矩阵论专用函数测试
# @File    : run.py
# @Software: PyCharm
from MatrixUtils import MatrixUtils
import numpy as np

if __name__ == '__main__':
    s = MatrixUtils(trace=True)
    a = np.array([
        [-1, 0, 1, 2],
        [1, 2, -1, 1],
        [2, 2, -2, -1]
    ])
    c = np.matmul(a.T, a)
    b = np.array([
        [5, 2, -4],
        [2, 1, -2],
        [-4, -2, 5]
    ])
    # print(L1)
    b = s.Hermite(a)
    print(b)
    # print(G)
    # print(c)
    # pass
